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Abstract 



The kinetic granular temperatures of a binary granular mixture in sim- 
ple shear flow are determined from the Boltzmann kinetic theory by using a 
Sonine polynomial expansion. The results show that the temperature ratio 
is clearly different from unity (as may be expected since the system is out 
of equilibrium) and strongly depends on the restitution coefficients as well 
as on the parameters of the mixture. The approximate analytical calcula- 
tions are compared with those obtained from Monte Carlo simulations of the 
Boltzmann equation showing an excellent agreement over the range of param- 
eters investigated. Finally, the influence of the temperature differences on the 
rheological properties is also discussed. 

PACS number(s): 45.70.Mg, 05.20.Dd, 51.10.+y, 47.50.+d 
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Typeset using REVTgX 



Although experimental and theoretical studies on granular media have been mainly fo- 
cused on assemblies of identical particles, there appears to be a recent growing interest among 
both theorists and experimentalists in the more complicated case of polydisperse systems. In 
this case, several kinetic theory studies in the freely cooling and thermostatted steady 
states have shown that the kinetic temperatures of each species are different. This violation 
of energy equipartition has been subsequently confirmed in experiments of vibrated granular 
mixtures [Q. When the system is sheared, a similar result has been recently found by Clel- 
land and Hrenya from molecular dynamics (MD) simulations of a binary-sized mixture 
of inelastic, smooth hard disks engaged in rapid shear flow. Their results were compared 
with previous kinetic theory calculations based on the assumption of equipartition of 
granular energy. As Clelland and Hrenya conclude , although this equipartition-of-energy 
assumption does not appear to have a negative impact on the ability of those earliest theories 
to predict the stress tensor in simple shear flow, a multi-temperature theory must be more 
appropriate. In the context of kinetic theory, the only primary attempts to include temper- 
ature differences in dense granular mixtures were put forward by Jenkins and Mancini 
and more recently by Huilin et al. [0. However, both works are phenomenological with no 
attempt to solve the kinetic equation. Instead, they assume that the velocity distribution 
functions are local Maxwellians. This can be reasonable to estimate the dense gas collisional 
transfer contributions to the fluxes, but not to compute their kinetic contributions. On the 
other hand, both theories are applicable to a general flow field. 

The goal of this brief report is to exphcitly get the dependence of the temperature ratio 
7 = T1/T2 on particle properties as well as on compositional parameters of a granular 
binary mixture subjected to the simple shear flow. The calculation of 7 allows one to assess 
the magnitude of the equipartition violation and its dependence on the parameters of the 
system. On the other hand, due to the complexity of the general problem, here we will 
restrict ourselves to the low-density regime in which case the velocity distribution functions 
fi for the two species verify a set of two coupled nonlinear Boltzmann equations. This set 
of equations is first analytically solved by using a first-Sonine polynomial approximation. 
Then, to check the reliability of our theoretical results, the direct simulation Monte Carlo 
method ^ is also employed to numerically solve the Boltzmann equation in the shear flow 
state. As will be shown later, our theory presents a much better agreement with simulations 
than the one given by Jenkins and Mancini 0. 

Let us consider a granular binary mixture composed by smooth inelastic disks or spheres 
of masses mi and m2 and diameters cxi and CT2. Collisions between particles are inelastic and 
characterized by three constant (independent) restitution coefficients an, 0:22, and ai2 = 021, 
where aij < 1 refers to the restitution coefficient for collisions between particles of species i 
and j. The mixture is under simple shear flow, namely, a macroscopic state with a constant 
linear velocity proflle U = a ■ r, where ake = aSkx^ey, a being the constant shear rate. In 
addition, the partial densities ni and the (global) granular temperature T are uniform. The 
time evolution of the temperature arises from the balance of two competing effects: viscous 
heating and collisional cooling. When both mechanisms cancel each other, the system reaches 
a steady state and the temperature achieves a constant value. This steady state is what 
we want to analyze here. From a microscopic point of view, the simple shear flow problem 
becomes spatially uniform in the local Lagrangian frame moving with the flow velocity U. 
In this frame, /j(r,v) — > /i(V), where V = v — U is the peculiar velocity. Under these 
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conditions, the set of Boltzmann kinetic equations read 

-aVy^MV) = J2J^AW^,a (1) 

X j 

where the Boltzmann colhsion operator Jij [V\fi, fj] describing the scattering of pairs of 
particles is 

[Vil/., /,] = af-' J dY, j da Q{a . gi2)(^ ■ gi2) [«-V.(v;)/,(v^) - /.(V0/,(V2^ 

(2) 

Here, d is the dimensionality of the system, aij = ((Xj + (jj) /2, is a unit vector along their 
line of centers, B is the Heaviside step function and gi2 = Vi — V2. In addition, the primes 
on the velocities denote the initial values {V^, V2} that lead to {Vi, V2} following a binary 
collision: V'^ = Vi - fiji (l + a^j^^j (ct ■ gi2)o- and = V2 + /ijj (l + a^^-^) (o- ■ gi2)o-. Here, 

fj,ij = mi/ {mi + mj). 

Our study is focused on the evaluation of the partial temperatures Tj, which measure 
the mean kinetic energy of each species. In terms of /j, they are defined as 



T, = J dY\m,V^U (3) 



The temperature of the mixture is T = XiTi + X2T2, where Xi = ni/{ni + 722) is the mole 
fraction of species i. The balance equation of the granular temperature for species i can be 
obtained by multiplying the Boltzmann equations by m^V^^ and integrating over V. The 
result is 

aPi^xy + ^CiPi = 0, (4) 

where pi = UiTi, 

P, = m,J dVVV/,(V), (5) 

is the partial pressure tensor of the species i and Q = J2j dj is the cooling rate for the partial 
temperature T,, with 

= / rfVm,yV.,[V|/„ /,] . (6) 

According to Eq. in the steady state the temperature ratio 7 = T1/T2 is given by the 
relation 

^ ^ X2(2Pl,xy ^ 
XlClP2,xy 

Thus, to get 7 one needs to determine the cooling rates Q and the xy element of the 
partial pressure tensors Pj. An equation for the elements of Pj follows immediately from the 
definition and the Boltzmann equation (|l]): 
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Ml 



(8) 



where we have introduced the colhsional moments A,, as 



A, 



mi 



dYVVJ,,[V\Uf, 



(9) 



The determination of Pj is a closed problem provided the moments Aij are explicitly known. 
This requires the knowledge of the velocity distribution functions /», which is quite an 
intricate problem, even in the elastic case. A useful way to estimate Q and Aij is to expand 
fi in Sonine polynomials. This approach is similar to the usual moment method for solving 
kinetic equations in the elastic case where the leading order truncation is known to be a good 
approximation. In the case of shear flow, we take the leading Sonine approximation, /i(V) — >■ 
/,,m(V) [1 + Q : D,(V)/2T,], where Q = (P,M) - U and D,(V) = m, [VV - iVyd)TL]. 
Here, 1 is the d x d unit tensor and f^M is a Maxwellian distribution at the temperature of 
the species i, i.e.. 



exp 



(10) 



With this approximation, the integrals appearing in the expressions of the cooling rates Cij 
and the colhsional moments A^ can be explicitly evaluated. Retaining only linear terms in 
Cj and after a lengthy calculation, one gets 



dT{d/2) 




1 + a 



2 - 



a,; 



9i + 9j 



6, 



3/2 



dT{d/2) 
9£ 



'1 + Oiij, 



X 



d 



d + 3 



+2- 



+ 2{d + 2) 9, " V 2(d + 2) 9, j ■ 



(12) 



In these expressions, Vq = y2T{mi +m2)/mim2 is a thermal velocity defined in terms of 
the temperature of the mixture T, and 



9i 



l + xi(7-l) 



9, 



l + xi(7- 1) 



(13) 



(14) 



Equations (|Tl]) and ([T^) extend previous expressions obtained in the three dimensional 
case. The approximation (p!2D allows one to solve the set of equations (P) and express Pj 
as a function of 7, while the approximation (^) gives the cooling rates Q. When all these 
expressions are used in Eq. (^, one gets a c/osed equation for the temperature ratio 7, that 
can be solved numerically. 
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A full presentation of the results is not possible because of the complexity of the pa- 
rameter space: an, 022, ai2, Tni/m2, Xi, and o"i/(T2. As in Ref. @], we assume for the sake 
of illustration that the spheres or disks are made of the same material {aij = a) and have 
the same mass density, i.e., mi/m2 = (cri/cr2)'^. This reduces the parameter space to a, Xi, 
and o"i/o"2- As expected, our results show that in general the kinetic temperatures of the 
mixture are different (7 7^ 1). There are only two trivial exceptions: the elastic case {a = 1) 
and the case of mechanically equivalent particles (mi = 1712, (J\ = (72). Beyond these cases, 
the dependence of 7 on the parameters of the problem is quite intricate. As an illustration, 
we plot the temperature ratio T1/T2 versus the diameter ratio 0"i/(j2 for a^i = | (Fig. [^) and 
Xi = \ (Fig. 1^) for three different values of a: a = 0.95, a = 0.9, and a = 0.8. We have 
considered the two dimensional case {d = 2). Also, for comparison we show the kinetic- 
theory predictions of Jenkins and Mancini 0. It is apparent that an excellent agreement 
between Monte Carlo simulations (symbols) and our theory is found over the entire range of 
values of size and mass ratios considered. Although the solid fractions considered by Clel- 
land and Hrenya 0] prevent us from making a quantitative comparison between our theory 
and their simulations, we observe that the behavior of 7 in dilute systems is qualitatively 
similar to the one found in Ref. Thus, for instance, at a given value of a the granular 
energy of the larger particle (say for instance, species 1) increases relative to that of the 
smaller particle as the size ratio 0"i/cr2 increases. Both Monte Carlo simulation and theory 
show that the temperature ratio presents a strong dependence on the restitution coefficient. 
With respect to the influence of composition, comparison of Figs. |l] and ^ indicates that 
7 exhibits a very weak dependence on the mole fraction xi. This behavior has been also 
found in recent experiments carried out on binary vibrated granular gases. We also see 
that all the above trends are qualitatively reproduced by the theory of Jenkins and Mancini 
1^ (which is restricted to nearly elastic disks), although these trends are however strongly 
exaggerated. Thus, for instance, at a size ratio o'i/a2 = 3, in the case Xi = \ and a = 0.8, 
the discrepancy between our theory and the simulation is less than 1% while it is around of 
115% in the theory of Jenkins and Mancini 0. 

An interesting point is to assess the influence of the temperature differences on the rhe- 
ological properties of the mixture. Although the comparison carried out in Ref. between 
previous theories (based on a single temperature) and simulation shows a quahtative good 
agreement at the level of the shear stresses, one expects that the new contributions coming 
from the energy difference leads to an improvement over previous theoretical predictions. In 
Fig. H, we plot the dimensionless stresses —P^y and P*y versus the restitution coefficient in 
the case xi = \ and mi/m2 = (^1/^2)^/^ = 10. Here, P* = P/nT = (Pi + P2)/nr. Also 
shown in this Fig. ^ is the result which would be obtained if the differences in the partial 
temperatures were neglected [i.e., 9i = jjijl in Eqs. (|lT]) and ([T^) ]. In general, inclusion of the 
two-temperature effects represents a significant improvement of the theory, especially in the 
case of the shear stress P*y, which is the most relevant rheological property in a sheared flow 
problem. This justifles the use of a two-temperature description to capture the dependence 
of stresses on dissipation. 

In summary, we have obtained an approximate evaluation of the temperature ratio of a 
sheared granular mixture. The accuracy of this calculation has been also conflrmed by Monte 
Carlo simulation of the Boltzmann equation. As was also found in recent MD simulations 
H], our results show that the temperature ratio strongly depends on dissipation and the 
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mechanical parameters of the mixture (especially on the ratios of mass and size). In addition, 
the effect of temperature differences on rheology is important (especially in the case of the 
shear stress) and leads to an improvement of the theoretical results with respect to the 
predictions made from a single-temperature theory. 

V. G. acknowledges partial support from the Ministerio de Ciencia y Tecnologia (Spain) 
through Grant No. BFM2001-0718. 
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FIGURES 



FIG. 1. Plot of the temperature ratio T1/T2 as a function of the size ratio oijoi = {mi/ 1712)^^"^ 
for a two-dimensional system in the case xi = 1/3 and three different values of the restitution coef- 
ficient a: (a) a = 0.95, (b) a = 0.9, and (c) a = 0.8. The sohd hues are the theoretical predictions 
while the symbols refer to the Monte Carlo simulation results. The dashed line corresponds to the 
prediction given by the theory of Jenkins and Mancini Q in the case a = 0.8. 



FIG. 2. The same as in Fig. [T| but for xi = 1/2. 



FIG. 3. Plot of the reduced elements of the pressure tensor P*y = Pyy/nT and P*y = Pxy/nT 
as a function of the restitution coefficient a for a two-dimensional system in the case mi/m2 = 10 
and xi = 1/2. The solid lines are the theoretical predictions while the symbols refer to the Monte 
Carlo simulation results. The dashed lines correspond to the theoretical results by assuming the 
equality of the partial temperatures T1/T2 = 1. 
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